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A combined mass and particle identification fit is used to make the first observation of the decay 
B° -> OfK"^ and measure the branching fraction of S° DfK^ relative to S° — » D+tt . This 
analysis uses 1.2 fb~^ integrated luminosity of pp collisions at = 1.96 TeV collected with the CDF 
II detector at the Fermilab Tevatron collider. We observe a B° — > OfK^ signal with a statistical sig- 
nificance of 8.1fT and measure B(B° -> DtK'f)/B(B° -D+tt") = 0.097 ± 0.018(stat) ± 0.009(sys). 
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One of the remaining open questions in flavor physics is 
the precise value of the angle 7 = arg {-VudV*,,/VcdV*^) 
of the unitarity triangle. Current measurements use the 
interference between b ucs and b cus diagrams 
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in B- -> d'^*^_^K^*^- and ^ D^*^°K^*^- decays 
when D*^ and -D'' are observed in common final states 
P, [3, H, 0, H, @], but suffer from the large difference 
between the amplitudes of these decays. With a large 
sample of hadronic decays, it may be possible to de- 
termine 7 from the interference, through _B"-_B^ mixing, 
of the same diagrams in the decay modes DfK~ 
and B2 ^ D-K+ 0, 1, which are expected to have 
a more favorable amplitude ratio; the two decays pro- 
ceed through color-allowed tree amplitudes whose ratio 
is suppressed by only a factor ~ 0.4 [§]. To determine 
7, a time-dependent measurement of the decay rates 
of 5° ^ D+K-, fiO ^ D-K+, 5° ^ DJK+, and 
DfK is required. The first steps in this effort 
are to observe these decay modes (which we will collec- 
tively refer to as B^ — > DfK^) and to measure the CP- 
averaged branching ratio B(B° DfK^) = \ [B{B° 
D+K-)+B{W, ^ D~K+)+B{B°, ^ D-K+) + B{B^, ^ 
D+K^)]. In this Letter we report the first observation 
of the B1 — + DfK^ decay modes and the first measure- 
ment of 6(5° -> DfK^)/B{W^ D+n-). We measure 
this branching fraction ratio since many of the systematic 
uncertainties cancel in the ratio and 0(5° Dfir^) is 
precisely measured elsewhere [Hi • 

We analyze pp collisions at = 1.96 TeV recorded 
by the CDF II detector at the Fermilab Tevatron collider 
with an integrated luminosity of 1.2 fb~^. A detailed 
description of the detector can be found elsewhere p^ . 
This analysis uses charged particle tracks reconstructed 
in the pseudorapidity [13| range |?7| < 1 from hits in a sili- 
con microstrip vertex detector [lj| and a cylindrical drift 
chamber [isj immersed in a 1.4 T axial magnetic field. 
The specific ionization energy loss [dE/dx) of charged 
particles in the drift chamber is used for particle identifi- 
cation (PID). A sample rich in bottom hadrons is selected 
by triggering on events that have at least two tracks, 
each with > 2 GeV/c and large impact parameter; 
the trigger further requires that these tracks originate 
from a secondary vertex well displaced from the primary 
interaction point . 

We reconstruct -6° Dfh^ candidates (where h — 
TT or K) as follows. First, we identify (l>{-^ 
K~ K~^)'K^ candidates [l3| using the invariant mass re- 
quirements 1013 < m{K-K+) < 1028 MeV/c^ and 
1948.3 < m{K-K+TT+) < 1988.3 MeV/c^. The D+ de- 
cay tracks must satisfy a three-dimensional vertex fit. 
No PID requirements are made on the Df decay tracks. 
We then pair the Df candidates with tracks to de- 
fine the B^ Dfh'f candidate sample, and require the 
Df-h~ pair to satisfy a three-dimensional fit to the B^ 
decay vertex. No mass constraint is used either on the 
(j) or on the candidate. Finally, we define a mass 
variable m(Z?s7r) for the DsTt hypothesis (i.e., assigning 
the daughter track h as a pion); m(Dsn) is used to pro- 
vide kinematic separation between the 5° DfK^ and 
DfTT~ signals. 

Further selection requirements are made to reduce 
combinatorial background. The discriminating variables 



are the transverse {\do\ < 60 pm) and longitudinal 
(jzo/czol < 3, where a^g is the uncertainty on zq) impact 
parameter of the B^ candidate with respect to the pri- 
mary event vertex; the transverse momentum (px > 5.5 
GeV/c) of the B'^ candidate; the isolation of the B'^ can- 
didate 



PT + E PT(track) ■ ' 

tracks 

w here the s um runs over tracks within AR — 
^/A(f>^ + Arf < 1 around the B^ direction originat- 
ing from the same primary vertex; the opening angle 
[AR{Df,h~) < 1.5] between the Df candidate and the 
track originating from the B'^ decay vertex (the "B^ 
daughter track"); and the projection of the and _D+ 
decay length along the transverse momentum of the 
candidate [L,y (5°) > 300 pm, L,y (fiO) /aL^^ (BO) > 8 
(where CFL^y is the uncertainty on L^y), and Lxy{Df) > 
Lxy (B°)]- The dE/dx calibrations are based on large 
samples of — > K~Tr'^ decays taken with the displaced 
track trigger. To avoid bias, the 5° daughter tracks are 
required to pass the same pt > 2 GeV/c trigger require- 
ment as the I?" — > K^n^ calibration tracks. 

Monte Carlo simulation is used to model signal and 
background and to determine trigger and reconstruction 
efficiencies. We generate single hadrons with BGEN- 
ERATOR [H, nil and simulate their decays with evtgen 
(20| . A detailed detector and trigger simulation is then 
performed. 

The greatest challenge in this analysis is to disentangle 
the various components contributing to the data sample. 
Apart from the S° DfK^ and D+n' signals, 

the sample contains partially reconstructed B^ decays, 
reflections from decays of other bottom hadron species, 
and combinatorial background. To separate the compo- 
nents and determine the number of candidates of each 
component type, we perform a maximum-likelihood fit. 
The fit variables are the invariant mass m{DsTr) of the 
candidate in the DsTt mass hypothesis and the PID vari- 
able Z, which is the logarithm of the ratio between the 
measured dE/dx and the expected dE/dx for a pion with 
the momentum of the B^ daughter track. The likelihood 
function is 

N M 

Hfi, ■ ■ ■ , fM-i) ^Yl^fjPj{m^)qj{Z,), (1) 

i=ij=i 

where = 1 — X]j=i^/j- The index i runs over the 
N candidates, and j runs over the M components; fj is 
the fraction of candidates in the jth component, to be 
determined by the fit. 

We group candidates into three categories by 
source. Combinations where the candidate and the 
track come from a single bottom hadron (_B°, B~ , 
A°) are called single-i? contributions. Non-bottom con- 
tributions where the candidate does not come from 
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a real Df are called fake--D+ combinatorial background. 
Combinations of a real Df with a track coming from 
fragmentation, the underlying event, or the other bot- 
tom hadron in the event are called real- 13+ combinatorial 
background. 

Mass probability density functions (pdf's) Pj{m) for 
the single- components are extracted from large sim- 
ulated samples of DfX decays, where the Df is 
forced to decay to 07r+. Separate mass templates are 
extracted for B° D+tt' and B° DfK^ fully re- 
constructed decays and for the partially reconstructed 
modes that overlap in mass with the — > DfK^: 
B° D+p~ and B° -> D*+tt-. Partially recon- 
structed modes missing more than one decay product 
are collected in one template. Contributions from the 
W D+{K-Tr+Tr+)X, -> D+ {K-Tr+Tr+)X and 
A|J ^ A+(pi^^7r+)X reflections are included. Likewise, 

we include Di*''^h^ decays (where h = Tr,K) 

whose branching fractions are known [2lj; the yields 
of these B'^ modes relative to each other are fixed to 
the values reported in [2ll |. Rather than parameteriz- 
ing the mass shapes, which are complicated for most 
of the single- i? components, we use histograms as mass 
pdf's. Sufficiently large Monte Carlo samples [approxi- 
mately 50 000 candidates after cuts of B^ Dfir^ and 
B^ — > DfK^] are generated to make the statistical fluc- 
tuations in the pdf's small. 

Special care has to be taken in the treatment of the low- 
mass tail of the decay mode B^ Dfn~{n-f), which is 
dominated by the radiative tail, and which overlaps with 
the DfK^ mass pdf. Improper accounting of the 

tail can bias both the measurement of the B^ -D+tt" 
yield and the B^ DfK^ yield by misidentifying a 
fraction of the B^ DfTT~ contribution as part of the 
(much smaller) B^ DfK^ contribution. The radia- 
tive tail is modeled in evtgen by using the PHOTOS al- 
gorithm for radiative corrections jl^] with a cut-off for 
photon emission at 10 MeV. We allow the normaliza- 
tion to float in the fit to account for uncertainties in 
the PHOTOS prediction of the size of the radiative tail. 
(The radiative tail of 5° — *■ DfK^ does not require 
special treatment. The kaon radiates less, and any re- 
sulting misidcntified B^ DfK^ contribution is easily 
absorbed by the other fit components, which dominate 
at m{D+TT-) below the DfK^ peak.) 

The mass distribution of the fake-D^ background is 
parameterized with a function of the form pbg(w) oc 
exp(— am) + /?. The shape parameters a and /3 are deter- 
mined in an ancillary mass-only fit of iJ^ candidates pop- 
ulating the sidebands of the D'^ mass distribution. To 
model the veal-Df background, we use a sample of same- 
sign DfiT^ candidates. A fit analogous to the fake-Z3+ 
case is performed on the D^7r+ mass distribution. Given 
to the limited statistics of the signal sample, we cannot 
separately resolve the real-D^ and fake-D^ combinato- 
rial backgrounds; in the default fit we therefore combine 
the two types of background. We assess a systematic 
uncertainty by allowing the relative size of the two back- 



ground types to vary. 

We determine the Z pdf's qj{Z) for pions and kaons 
from D*+ D"{K-TT+)7r+ decays. The flavor of 
the daughter tracks of the D'^ in the decay D*+ 
£)°(i<f^7r+)7r+ is tagged by the charge of the soft pion 
from the decay. Taken together with the large 

signal-to-background ratio of the Am = m{K~TT~^Tr'^) — 
m[K~TT^) peak, this charge tagging yields a very clean 
sample of pions and kaons. We further reduce back- 
ground contamination by sideband-subtracting in Am. 
The mean values of Z for kaons and pions are separated 
by approximately 1.4 standard deviations. The Z distri- 
butions for both species (shown in Figure [T]) have sim- 
ilar widths. Because the data sample contains semilep- 
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Figure 1: Z distributions for pions and kaons derived from 
prompt Z)*+ decays. 

tonic decays, we need to model the dE/dx distributions 
of muons and electrons as well. For muons, which are 
a small contribution in the mass region of interest, the 
pion template can be used without introducing a signif- 
icant systematic uncertainty. For electrons, we derive a 
template from a J/ip — > e"^e~ sample. The Z pdf for the 
fake-D^ combinatorial background is determined from 
data by selecting candidates from the sidebands of the 
Df mass distribution. All Z pdf's are represented as 
histograms. 

Figures [5] and [3] show the fit projections in mass and 
Z. The yields determined b_y the fit are 1125 ± 87 
D+TT- and 102 ± J_8 S° DfK^ candidates. 
The branching fraction of B^ — > DfK^ relative to B^ 
D^TT^ , corrected for the relative reconstruction efficiency 
e^/eK = 1.071 ± 0.028, is ^(5° -> DfK^)/B(B°, 
D+TT-) = 0.097 ± 0.018. A fit performed with the 
B° DfK^ yield set to zero is worse than the default 
fit by AlogL = —32.52; the corresponding statistical 
significance of the — > DfK^ signal is 8.1 standard 
deviations. 

Systematic uncertainties on B(B° DfK^)/B{Bl 
DfTT~) are studied by incorporating each effect in the 
generation of simulated experiments which are then fit- 
ted using the default configuration. The bias on B{B'^ 
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Table I: Systematic uncertainties on B{B^l 
DfK^)/B{Bl-^Dt^-). 



Source 


Systematic uncertainty 


dE/dx pdf modeling 


0.007 


Mass pdf modeling 


0.004 


Combinatorial-background model 


0.002 


Fitter bias due to finite statistics 


O.OOf 


Sum in quadrature 


0.009 



DfK'^)/B{B° D+TT-), averaged over 10000 simulated 
experiments, is used as the systematic uncertainty asso- 
ciated with the effect under study. Table [T] summarizes 
the systematic uncertainties. The systematic uncertain- 
ties are dominated by the modehng of the dE/dx (0.007), 
specifically by the differences between the Z distributions 
of daughter tracks (from which the kaon and pion Z 
pdf's are derived) and Bl daughter tracks; these differ- 
ences arise from effects such as the greater particle abun- 
dance in the vicinity of a prompt D*^ compared to a iJ^, 
and hence a higher probability for daughter tracks to 
contain extraneous hits. Modeling of the mass distribu- 
tions of the single- i? components (0.004), which includes 
statistical fluctuations in the mass pdf's, and modeling 
of the combinatorial-background mass shape (0.002) are 
comparatively minor contributions. The total systematic 
uncertainty is obtained by adding the individual system- 
atic uncertainties in quadrature; at 0.009, it is about half 
as large as the statistical uncertainty. 

The analysis procedure was crosschecked in several 
ways. Most importantly, before performing the mea- 
surement on the — > DfK^ signal sample, we veri- 
fied our method using two control samples, i?" — ^ D~^X 
and D*~^X. In both cases, our results are statis- 

tically consistent with world-average values. We mea- 
sure B{W -> D+K-)/B{W D+TT-) = 0.086 ± 
0.005(stat), 1.0 standard deviations from the world aver- 
age; and B(B^ D*+K-)/B{W ^ D*+t:-) = 0.080 ± 
0.008(stat), 0.3 standard deviations from the world av- 
erage 23]. The relative branching fractions B{B'^ 
D+p-)/BiW ^ D+n-), BiW ^ D*+TT-)/B{W ^ 



D+TT-), and B{B° 



D*+p-)/B{B^ 



D 



') were 



also found to be consistent with world averages. Finally, 
the fractional size of the radiative tails of i?'^ — > D+tt", 



B° 



D*+TT-. and B2 



Dtir are found to be in 



agreement with each other (and about twice as large as 
the PHOTOS prediction). 

In conclusion, we have presented the first observation 
of the -B^ DfK^ decay mode with a statistical sig- 
nificance of 8.1 standard deviations. The -8° DfK^ 
event yield is 102 ± 18 (statistical uncertainty only). We 
use this sample to measure B{B^ DfK^)/B{B^ 
D+TT") = 0.097 ±0.018(stat)±0.009(sys). This resuh is 
consistent with naive expectations based on the branch- 
ing fraction ratio for the analogous and B~ decays, 
taking into account also the expected contribution from 
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Figure 2; Mass projection of the likelihood fit. Fit com- 
ponents are stacked. B DX denotes all single-B con- 
^ Dt-K' and ~* DfK^, namely 
D*'^-R~, partially reconstructed 
missing more than one decay prod- 
+ 7r+)X, B" ^ D+{K--K+7v+)X and 



tributions except 



B's 

net 
A? 



D, 



X modes 
• D+{K--n 



Kt{pK~-K^)X reflections, and B° 



dI*'' K'^; the small peak in the B 



and 

DX template 

is due to the {K~ tv'^tt^)tv~ reflection. The resid- 

ual plot at the bottom shows the discrepancy (data minus 
fit) in units of standard deviation (a); for the bins with low 
statistics, neighboring bins are combined until the predicted 
number of events is greater than five. The of the projection 
is 79.0 for 72 degrees of freedom. 



^ DjK+ decays ^. 
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Figure 3: Z projection of the likelihood fit in the region of 
interest for B° ^ D^K^ (5.26 < m{Ds-n:) < 5.35 GeV/c^). 
Fit components are stacked. The residual plot at the bottom 
shows the discrepancy (data minus fit) in units of standard 
deviation (a); for the bins with low statistics, neighboring bins 
are combined until the predicted number of events is greater 
than five. The of the projection is 30.7 for 14 degrees of 
freedom. 
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